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Abstract. 

In this paper we discuss the demographics of the radio blazar popu- 
lation: (i) what are their parent ('unbeamed') sources and (ii) what mag- 
nitude and/or type of evolution have they undergone ? The discussion is 
based on models of radio source evolution and beaming based on a 'dual 
population' unification paradigm. These models, developed from radio 
blazar properties in bright samples, predict blazar demographic trends at 
the lower flux-density levels; samples from deep mJy- level surveys {e.g. 
NVSS and FIRST) may now provide direct tests of these predictions. 



1. Dual-population Unification 

All galaxies are radio galaxies in a sense. The scale of activity ranges from 
'normal' galaxies (like our own) whose radio emission emanates from supernova 
remnants. At higher radio powers are the 'starburst' galaxies, greatly-enhanced 
stellar activity resulting in stronger radio emission. At similar radio powers are 
the Seyfert galaxies with their small-scale radio core-jet structures and with clear 
indication of AGN activity. At the top end of the scale are the powerful radio 
galaxies, together with BLLacs and quasars, the latter two being collectively 
termed blazars. 

Radio galaxies tend to have steep radio spectra and therefore dominate 
low radio-frequency (y < 0.5 GHz) surveys. The structures are predominantly 
double-lobed; Fanaroff & Riley (1974) compared the distances between cen- 
tral maxima of the radio structures to overall size and discovered a dichotomy. 
Sources with regions of maximum brightness separated by more than 0.5 times 
the overall source size are edge-brightened (FRII), while when the regions of 
maximum brightness lies within this limit the source has diffuse outer lobes 
(FRI) (centrally concentrated); and the highest radio power objects are predom- 
inantly FRIIs whilst those of lower radio power are usually FRIs. The overlap 
between the classes at middling radio powers is seen in the local radio luminos- 
ity function (Figure 1), and the division is known to be a function of optical 
luminosity as well as radio (Owen and Ledlow 1993). 
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Figure 1. The local radio luminosity function at 151 MHz (dark, solid 
line) comprises contributions from 3 populations: starburst galaxies 
(dotted), FRI (dot-dash) and FRII radio galaxies (dashed). 



Evidence is accumulating that the powerful radio sources, double-lobed ra- 
dio galaxies on the one hand and blazars on the other, are 'unified' by projection 
effects, with a central torus together with relativistic beaming responsible for 
our classification of objects as radio galaxies or blazars (Scheuer & Readhead 
1979; Orr k Browne 1982; Scheuer 1987, Barthel 1989, Urry & Padovani 1995 
and references therein). The 'dual-population' unified scheme (Wall & Jackson 
1997), summarised in Table 1, posits that FRII radio galaxies are the parent pop- 
ulation of all quasars and some BLLac type objects. The FRII-quasar sources 
are galaxies with class 'A' spectra (Hine & Longair, 1979), having strong, high- 
excitation emission lines in their optical/UV spectra. The FRII-BLLac sources 
are galaxies with class 'B' spectra (Hine & Longair, 1979), having only weak, if 
any, high-excitation lines in their spectra. The second part of the scheme posits 
that FRI radio galaxies are the parent population of the remainder of the BL Lac 
type objects. 



2. The cosmic evolution of the parent populations 

Our analysis interprets radio source counts and and identification/redshift data 
for samples at high flux densities, adopting a dual-population unified scheme 
paradigm. Using a simple parametric model to describe the evolution of the 
underlying parent populations we first fit low-frequency radio data (v < 0.5 
GHz), where radio samples are unbiased by the effects of Doppler beaming. 
Then, using these evolution models for the parent populations of FRI and FRII 
objects, we fit the 5 GHz radio source count using a small number of parameters 
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Table 1. Unified classes of powerful extragalactic radio sources 



UV/optical Radio 
emission spectrum Source 
Population type at 5 GHz type 



FRII, class A 



narrow 
broad 
broad 



steep 
steep 
'flat' 



RG 
RG 



Quasar 



FRII, class B 



weak/none 
none 



steep 
'flat' 



RG 
BL Lac 



FRI 



weak/none 
none 



steep 
'flat' 



RG 



to describe the Doppler beaming, which gives rise to blazars (i.e. quasars and 
BL Lacs with c^tgHz > -°- 5 > where S « 

The fundamental dataset to our analysis is the differential radio source 
count, the count at a single observing frequency of source density on the sky as 
a function of flux density. Compiling data from wide-area and deep pencil-beam 
surveys yields source counts spanning 2 decades in frequency and 6 decades in 
flux density as shown in Figure 2. 

In fact radio source counts contain a wealth of cosmological information, 
in particular providing potent evidence against a steady-state Universe (Ryle & 
Clarke 1961) and indicating that differential evolution must take place, with the 
most luminous sources showing the strongest evolution. In addition the shape 
of the source count changes with observing frequency, suggesting that more 
complex behaviour than just number density evolution is taking place - and our 
analysis shows just what this behaviour is. 

Identification programmes and further studies of radio source samples have 
revealed distinct regions within the counts. These regions are indicated along 
the abscissa in Figure 2 and discussed briefly here: 

Region 1 At the very highest flux densities the source count is near-Euclidean, 
although typically less than 20 sources contribute to this region. This 
region is comprised of a mixture of nearby sources and cosmologically- 
distant powerful sources. The flat count arises due to the dilution of the 
powerful evolving sources by the local lower-radio-power sources which are 
far more abundant at low redshifts. 

Region 2 In this region the counts rise more rapidly with decreasing flux 
density than the Euclidean (-3/2 power law) prediction, and then reach a 



1 Note that our definition of blazar is solely a radio definition and says nothing about the optical 
equivalent width criteria. 
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Figure 2. Source counts in relative differential form AN/ANq where 
AN is the number of sources per sterad with flux density S u between 
S2 and Si and ANq is the number of sources expected in a uniformly- 
filled Euclidean universe (N = K u Sv 3/2 with K u = 2400, 2730, 3618, 
4247, 5677 and 3738 for the six frequencies shown). The horizontal 
bars show the flux-density bin width, S2 to Si and the vertical error 
is the s/N error. The curves are polynomial least-square fits to the 
counts. The counts are compiled from radio survey data described in 
Wall (1994). 
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Euclidean plateau. This region of the count is made up of the most power- 
ful radio sources at cosmological distances. The low-frequency counts (< 
0.5 GHz) comprise almost entirely of steep-spectrum sources whose radio 
emission is greatest at these frequencies. The width of the plateau increases 
with increasing frequency, due to an increasing blazar (flat-spectrum) con- 
tribution. 

Region 3 Below the Euclidean plateau the counts at all wavelengths drop 
away from the Euclidean prediction. This region extends up to two or- 
ders of magnitude in flux density. (Limited) identification data shows that 
the sources at such flux densities are lower-power objects at intermedi- 
ate redshifts and not, as might be supposed, powerful radio sources at 
ever-increasing redshifts. Very little work has been done on the blazar 
population at these flux density levels. 

Region 4 At low flux densities the counts again flatten to near-Euclidean. In 
this region the dominant populations change dramatically to 'blue' star- 
burst galaxies and 'red' FRI-type galaxies, both relatively local and of low 
radio powers (e.g. Windhorst et al. 1993, Benn et al. 1993). 

Using the spectroscopically-complete 3CRR sample at 178 MHz (Laing, 
Riley & Longair 1983 and more recently published data collated by R. Laing, 
private communication) we first determine an appropriate form for the evolution 
of the FRI and FRII populations using the (V/V max ) statistic. 

Figure 3 shows V/V max values for 137 steep-spectrum FRII radio galax- 
ies and 26 FRI radio galaxies. Both populations show increasing (V/V max ) 
with radio power. However, FRIs exibit negative evolution with {V/V max ) 
< 0.5 at log 10 (Pi5i MHz ~ 25.0) and the FRIIs exhibit strong positive evo- 
lution with (V/V max ) > 0.7 at log 10 (Pisi mh z ~ 27.5). Interestingly, there is a 
trend for (V/V max ) to continuously increase with log 10 (Pi5i mh z ) across the pop- 
ulations, suggesting luminosity-dependent evolution which may be population- 
independent. 

Adopting exponential luminosity-dependent density evolution (LDDE) we 
fit the source count at 151 MHz comprising the 3CRR and 6C survey (Hales, 
Baldwin, & Warner 1988), using a minimisation routine. The model fit requires 
strong evolution of the most powerful FRII sources, with little or no evolution 
of the FRI population (Table 2). Thus the fit reproduces the {V/V max ) be- 
haviour for the radio galaxy population as a whole, in that only FRIIs with 
log 10 (.Pi5i MHz) > 25.44 undergo any evolution, the lower-power FRIIs and all 
FRIs have a constant space density to their cut-off redshift, z c . 

3. Blazars from beamed parent radio sources 

We fit the 5 GHz source count, which is well-defined over a wide flux-density 
range, starting from the evolution model fit of Table 2. However at 5 GHz we 
incorporate the beamed products (blazars) by randomly aligning the sources 
with respect to our line-of-sight, allowing both different Lorentz factors (7) and 
different core-to-extended flux ratios for the two (FRI, FRII) populations. The 
best-fit beaming parameters are shown in Table 3. 
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Figure 3. V/l/ma^ for individual 3CRR FRI radio galaxies (*) and 
FRII radio galaxies (o). Overplotted are (V/V max ) values for bins of 
0.5 in log 10 (Pi5i MHz) for FRIs (dashed +) and FRIIs (solid +). 



Table 2. Fitted evolution parameter values at 151 MHz 
for LDDE: p(P,z) = p {P) exp M(P)t{z) 

Exponential LDDE Chi-square test 

Population parameter values Xmin u<1 



FRI 


M max = 0.0, z c ■- 
Pi k, P2 not used given 


= 5.0 

M max = 0.0 




FRII, class A & B 


M max = 10.93, z c 
P l = 25.44, P 2 = 


= 5.62 
27.34 






best-fit 


30.73 


33 



a Degrees of freedom 
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Table 3. Fitted beaming parameter values at 5 GHz 

Chi-square test 

Population Parameter values Xmin ^ a 



FRI 



7 = 15.0, R med oc P, 



>-0.55 
151 MHz 



FRII, class A & B 



7 = 8.5, Rm ed = 0.01 
^c(Rmed) = 7°.l 



best-fit 



32.98 



25 



The contribution from each population to the 5-GHz source count is shown 
in Figure 4. The contribution from the quasar population peaks at a higher flux 
density limit than its parent population, and a similar behaviour is seen for the 
the BL Lac objects. 

We can analyse the relative importance of each class of radio source as 
a function of flux density. Blazars are dominated by quasars (high-excitation 
FRII parents) at S5 GHz > 0.3 Jy, but dominated by BL Lac-type objects at lower 
flux densities. The magnitude of the rapidly changing quasar fraction predicted 
by our model is well-matched by observations from the 0.25 Jy-sample of flat- 
spectrum sources of Shaver et al. (1996) - see Figure 15 of Jackson Sz Wall 



In terms of population mix, it is also instructive to consider the integral 
population count, shown in Figure 5 in which the model fit has been transposed 
to 1.4 GHz. This figure reveals that complete samples are predicted to comprise 
~20% BL Lac objects for flux density limits extending down to Si.4GHz ~ 0.1 
mJy Optical spectroscopy of objects selected from the FIRST survey (Becker, 
White & Helfand 1995) could provide substantial tests of this prediction to 
£1.4 GHz ~ 1 mJy. 

4. Redshift distributions of blazars 

Transposing our model (evolution plus beaming) to 2.7 GHz we compare our 
predictions with the limited observational data currently available. Figure 6 
shows the mean redshift, (z), across a wide flux density limit range for all source 
types as well as three sub-classes. For complete samples the model predicts that 
the value of (z) peaks at ~20 mJy then drops rapidly. In detail, we find that 
(i) the contribution from steep-spectrum sources dominates and therefore traces 
the overall (z), (ii) for the quasar population, (z) rises continuously as S2.7 GHz 
decreases. However the number of quasars is negligible below 52.7 GHz = 20 mJy 
so this has little effect on the overall (z) value and (Hi) the peak (z) of the BL 
Lac sources occurs at a higher flux density than that of the quasars, although 
the value of (z) is much lower. 

Comparing the model (z) values to the 0.25 Jy flat-spectrum quasar sample 
of Shaver et al. (1996 and in preparation) we find reasonable agreement. The 
data points for the 0.25 Jy sample shown on Figure 6 are lower limits as this 



(1999). 
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5-GHz Flux Density / Jy 



Figure 4. Observed (+) and model (solid line) differential source 
counts at 5 GHz. The model count comprises 4 contributions from 
FRII parents (dashed): 1) high-excitation radio galaxies, 2) quasars 
(from high-excitation parents), 3) low-excitation radio galaxies and 4) 
BL Lac type objects (from low-excitation parents). There are 2 con- 
tributions from FRI parents (dotted): 5) radio galaxies and 6) BL Lac 
type objects. The contribution from the starburst galaxy population 
is shown dot-dashed (7). 
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Figure 5. Integral population mix from the model fit at 5 GHz trans- 
posed to 1.4 GHz. Numbers in braces refer to the populations described 
in Figure 4. 
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Figure 6. Model (z) of all sources to limiting flux density at 1.4 GHz 
(solid line). Contributions from the steep spectrum FRII, FRI and 
starburst galaxies (populations 1, 3, 5, and 7) shown dotted, quasars 
(from high-excitation FRIIs, population 2) dot-dashed and BL Lacs 
(from low-excitation FRIIs and FRIs, populations 4 and 6) dashed. 
Data points are from the flat-spectrum quasar sample (*) (Shaver et 
al. 1996) and the 2 Jy sample (o) (Wall & Peacock 1985, Morganti et 
al. 1997). The solid curve shows the total. 

sample selects quasars with (a^^GHz ^ —0.4). Our model, however, defines 
blazars as having (cc^ghz > —0.5). Beamed objects with a^GHz > —0.4 have 
parent sources of lower radio luminosity, which in turn are predicted to have 
a lower (z) distribution than the less-beamed —0.5 < a^GHz < —0.4, higher- 
luminosity objects which are missing from the 0.25 Jy sample. 

The (z) value for the complete 2-Jy sample (Wall & Peacock 1985, Morganti 
et al. 1997) agrees well with the overall (z) model prediction. 

Figure 7 splits the BL Lac (z) distribution of Figure 6 into its contributing 
parent populations. We see that low-excitation-FRII BL Lac sources have a (z) 
distribution which is very similar to the quasar distribution. However, the FRI 
BLLacs have a very different distribution, with (z) gradually increasing with 
decreasing flux density. Morphological studies of a sizeable sample of BLLac 
objects from e.g. the FIRST survey would provide direct, powerful tests of this 
prediction. 
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